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Квадратичные оптимизационные задачи 
компьютерной геометрии 


Работа посвящена постановке и решению класса квадратичных оптимизационных задач компьютерной 
геометрии: поиск эллипсоида минимального объема, содержащего множество точек евклидового пространства, 
поиск минимального расстояния между эллипсоидами, построение гиперплоскости, разделяющей два 
эллипсоида. Предложены эффективные алгоритмы для решения этого класса задач. 


Введение 


Компьютерная, или вычислительная, геометрия - это новое направление иссле- 
дований, предметом которой является анализ и построение эффективных алгоритмов 
решения геометрических задач [1]. Она рассматривает такие актуальные задачи, как: 
триангуляцию, построение выпуклых оболочек, определение принадлежности одного 
объекта другому, поиск их пересечения, разделение множеств, аппроксимацию и т.п. 
Компьютерная геометрия используется в распознавании образов, инженерном проекти- 
ровании, медицине, компьютерной графике, робототехнике, компьютерных играх и т.д. 
Большинство рассматриваемых задач в этих областях являются оптимизационными, 
для решения которых разработано большое число методов. Однако решение таких за- 
дач в режиме реального времени требует разработки новых и эффективных алгоритмов. 

В последние годы внимание исследователей в этой области привлекают такие 
объекты, как эллипсоиды. Эллипсоиды имеют небольшое количество геометрических 
параметров и эффективны для приближения широкого класса выпуклых объектов 
при моделировании физических систем. Обнаружение столкновения или наложения 
двух эллипсоидов является, таким образом, важной проблемой с приложениями в ком- 
пьютерной графике, компьютерной мультипликации, виртуальном мире, робототех- 
нике, САО/САМ, вычислительной физике. Однако использование эллипсоидов огра- 
ничивалось недостатком эффективных методов для обнаружения их разделения или 
наложения. Рассмотрение эллипсоидов в и-мерном евклидовом пространстве является 
стимулирующим фактором для разработки эффективных алгоритмов квадратичной 
оптимизации. 


Постановка задач и анализ последних результатов 


Любой геометрический объект можно представить множеством точек {х'} в п- 


мерном евклидовом пространстве. Эти точки можно описать эллипсоидом минималь- 
ного объема. Уравнение эллипсоида имеет вид 


Ехо) = {х|(«-х0)Г0(х-х0)< 1}, 
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где точка хо — центр эллипсоида, а О-— положительно определенная матрица. Известно, 


что эллипсоид минимального объема, содержащий заданное множество точек, един- 
ственный. Объем эллипсоида равен 


КЕ) = ——® 


\4еко) 


где @- объем единичного шара. Таким образом, задача построения эллипсоида, со- 


— оаеко-®)"?, 


держащего заданное множество точек 4х, сводится к решению следующей оптими- 
зационной задачи: 


ши {4еко") (хх) Ох -х0) < Ы=1,...т,О > 0}, (1) 
где необходимо определить положительно определенную матрицу О и вектор хо 


Учитывая то, что функция ае(О-') является выпуклой, а множество всех положитель- 


но определенных матриц образует выпуклый конус, задача (1) является задачей вы- 
пуклой оптимизации. Однако условие положительной определенности матрицы не мо- 
жет быть записано в явном виде, что создает трудности в решении задачи (1). Заменим 


целевую функцию задачи (1) на эквивалентную 1ю8(4еКО')). Градиент этой функции 
равен О. Это позволяет использовать для решения задачи (1) метод Франка-Вулфа, в 


котором на каждой итерации решается линейная задача полуопределенной оптими- 
зации 

ши{О, хО | (хх) О(х'- м) Ы=Ь,...т‚О > 0}. 
В настоящее время разработаны прямо-двойственные методы внутренней точки для 
решения линейных задач полуопределенной оптимизации [2]. Ниже рассматривают- 
ся более простые алгоритмы решения задачи (1). 

Если геометрический объект задан в виде многогранника {х | Ах < 0,х> 0}, то 
задача построения эллипсоида минимального объема, содержащего многогранник, 
значительно усложняется, так как становится МР-сложной [3]. 

Для поиска минимального расстояния между двумя эллипсоидами разработан 
алгоритм [4], в котором на каждой итерации строятся вписанные в эллипсоиды шары 
максимального радиуса, касающиеся граничной точки эллипсоида. Новое приближе- 
ние искомого отрезка находят на прямой, соединяющей центры шаров. Другие алго- 
ритмы на каждой итерации решают вспомогательную квадратичную задачу поиска 
минимального расстояния от точки х * до границы эллипсоида 


пах х* |7 (х-х0)ГО(х- 0) < В. (2) 


В работе предлагается эффективный алгоритм для решения задачи (2). 

Проблема классификации объектов — одна из наиболее важных в области искус- 
ственного интеллекта. Разработано множество алгоритмов для построения гиперплос- 
костей, разделяющих различные выпуклые объекты [5], [6]. Разделяющую гиперплоскость 
легко построить, если найдено минимальное расстояние (отрезок) между эллипсоида- 
ми. Она будет перпендикулярной минимальному отрезку и проходить через его середину. 


Алгоритмы решения задач компьютерной геометрии 


Рассмотрим задачу построения эллипсоида минимального объема, содержащего 
заданное множество точек. Очевидно, что точки, принадлежащие внутренности их 
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выпуклой оболочки, можно исключить при построении эллипсоида. Для исключения 
таких точек решим последовательность задач линейного программирования 
ми Уж. => 
ЛЕК) ЛеЦЮ ЛЕЦ®) 
для К =1...., (К), где 1(К)— множество индексов оставшихся точек. Если 0 < тах{Я,} <1, 


то точку х^ можно представить линейной комбинацией других точек, тогда такие точ- 


ки исключаем из дальнейшего рассмотрения (обновляем множество индексов [(К)). 
Пусть число оставшихся точек равно ит. 

Разработан простой алгоритм решения задачи (1), когда оси эллипсоида совпа- 
дают с осями координат [7]. В этом случае, задача (1) эквивалентна следующей задаче: 


п п 
тах {> `102(а,)[У`а;(х, -х0)* <Ы =. .тьа > 0}. (3) 
1 = 
Задача (3) преобразуется к двойственной 
п 
шах {> `102(и‚(2)-у7(2))[е'2=12>0}, (4) 
= 


где е - вектор, состоящий из единиц, а 


т 
ира) Е Мат, 
1=1 


т 
у,(2) = Ух] ый 
1=1 


Для решения задачи (4) используем метод Франка-Вулфа [8]. Пусть 2*- ее решение, 
тогда искомый эллипсоид имеет вид 


2 
м (х У (2 )) 
* 2(7*)) _ 
ит (и (=*) 5 (2*)) 
Для построения эллипсоида, в общем случае, преобразуем систему координат, выби- 
рая новые координатные гиперплоскости методом наименьших квадратов 


Е = {х 


т 
: в. 2 
ши {> (ах -Б) а 1, 
1 
т 
где а х=Ь- искомая гиперплоскость. Следующую координатную гиперплоскость 
найдем методом наименьших квадратов для проекции точек на пересечение преды- 
дущих найденных координатных гиперплоскостей. Проекция точки х’ на пересече- 
ние гиперплоскостей {х | Ах =Ь} вычисляется по следующей формуле: 
х=х! — АГ(ААТ) ЦАХ-Ь). 
После преобразования координат решаем задачу (4) и для найденного эллипсоида 
производим обратное преобразование координат. 


Разработана программа, реализующая данный алгоритм. Пример работы програм- 
мы представлен на рис. 1 для 7 точек 
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Уравнение эллипсоида имеет вид 


0,0103(х, +0,641)? +0,0367(х, 0,997)? +0,0064(х, +0,641)(х, -0,997) =1. 


г 


Рисунок 1 — Эллипсоид, содержащий точки 


Задаче поиска минимального расстояния между двумя выпуклыми множества- 
ми 51 и 5› посвящено много работ [9], [10]. Она формулируется следующим обра- 


зом. Найти 


п {| х-у||хезьуе52}. 
Будем предполагать, что заданы два эллипсоида матрицами О и О), с центрами 


в точках 21,22 соответственно. Для поиска минимального расстояния между эллип- 


соидами Е (2') и Е>(22 ) рассмотрим алгоритм, предложенный в [4]. 


1. Находим точки х',х” пересечения отрезка [2,22 ] с границами эллипсоидов 
Е,Е>. 
2. Вычисляем градиенты 81,2. в точках х',х? к эллинсоидам Е, Е›. Задача по- 
иска минимального расстояния решена, если 
РИ 5 Ге 
(2-х) (е-—т в) = =12, 
151 
| 
тогда ее решение отрезок [х ‚х’] (е- вектор, состоящий из единиц). 
Определяем центры новых шаров по формулам 
О 
й 


и переходим к шагу 1. 
Из сходимости данного алгоритма [4] следует, что скалярные произведения 


Г _* 
(2:1) 8; ; 

К * > 
|=; ||| 2; || 


* 
будут убывать, когда К —> ©, где ©; — градиент в оптимальной точке. Это означает, 


что скорость сходимости алгоритма можно увеличить, если на втором шаге выбирать 
К К 
градиент в виде ©; +1442; , где д, 0. 


В других алгоритмах на каждой итерации решается только задача (2). Она может 
быть решена методом Франка-Вулфа. Более простым является следующий алгоритм: 
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1. Находим градиент 2> (х') и решаем задачу 
тах {27 х|хеЕЕ.}\. (5) 


Ее решение достигается в точке 


\ 52 @> 8> 
2. Проектируем точку х! на гиперплоскость р) (х-—х*) =0, имеем 
Тож. 1 
х*-—х 
Е ) 
|182 || 


Если точки х и х* не совпадают, то полагаем х' = (х+х*)/2 и переходим к 


82. 


шаху 1. 
3. В противном случае, проектируем исходную точку на гиперплоскость 


= (х-х*) =0. Если проекция хР совпадает с точкой х*, то задача решена, иначе по- 
лагаем х' = (1-а)х'+ах?(0<а<Ъа — 0) и переходим к шагу 1. 


На каждой итерации этого алгоритма расстояние между точкой хи границей 
эллипсоида Е> убывает. Одновременно убывает расстояние между точкой касания 
гиперплоскости эллипсоида и проекцией на эту гиперплоскость исходной точки. Поэто- 
му в пределе будет найдено решение задачи (2). 

Были проведены многочисленные эксперименты по сравнению эффективности 
этих алгоритмов. Они показали, что скорость сходимости второго алгоритма выше. 
В примере (рис. 2.) 

Е ={х|1,2(х 2)? +0,05(х-3)* < В, 
2 2 
Е> == {х (2х +3) + 20х> + 4(х! +3)х> < 103, 


первый алгоритм потребовал более 300 итераций (нахождения шаров) для заданной 
точности &=0,00001, а второй алгоритм потребовал 7 решений задачи (2). 


Минимальный отрезок равен [1,348937,—0,13478;0,531722,-—0,32035]. 


=. 9 


Рисунок 2 — Минимальный отрезок, соединяющий эллипсоиды 
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Рассмотрим задачу построения гиперплоскости, разделяющей два эллипсоида. Для 
ее решения используем алгоритмы предыдущей задачи. Если [х',х2] — кратчайший 
отрезок между двумя эллипсоидами, то разделяющей гиперплоскостью будет следующая 

Ре 
х-х 


2 


Эта гиперплоскость перпендикулярна отрезку [х,х?] и проходит через его середину. 


(д -хоГ(х- )=0. 


Выводы 


Рассмотрены задачи компьютерной геометрии, которые относятся к квадратич- 
ной оптимизации. Предложены эффективные алгоритмы для решения этого класса задач. 
Их сходимость следует из сходимости используемых алгоритмов оптимизации. Разра- 
ботано программное обеспечение, реализующее данные алгоритмы, и проведены чис- 
ленные эксперименты. Эти эксперименты подтверждают эффективность предложенных 
алгоритмов. 
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А.1. Косолап 

Квадратичн! оптимзащйн! задач! комп’ютернот геометри 

Робота присвячена постановцщ! та розв’язку класу квадратичних оптимзащйних задач комп’ютернот 
геометри: пошук елшсойду мимального об’ему, що мстить множину точок евкллдового простору, 
пошук м!ин1мально] в1дстан1 мж елшсо1дами, побудова гшерплощини, що роздляе два елшсо1да. 
Запропонован! ефективн! алгоритми для розв’язку цього класу задач. 


А. Коздар 

Опадгайс Орйпи7аНоп Ргоет$ оГ Сошршег Сеотейгу 

ТБе рарег 15 4еухфе4 ю Фе зетепе ап фе зоавоп оЁа с1а5$ оЁ ападгайс орипияше ргоепт$ ш Фе 
сотрщег оеотебу: Фе зеагсВ оЁ бе тшипит уо[лте еШрзо!4 фаё сотшатз Фе зе! оЁ рош5 ог ЕисИаеап 
зрасе, фе зеагсВ оЁ фе пшипит 41апсе Бебуееп еШрзо14$, пе сопзгасНоп оЁ Фе Вуреграпе зерагайпе мо 
ерзо14$5. ТВе еНесйуе а[еогИилз$ Гог Пе зоШоп оЁ 1$ с1а$$ оРргоет$ аге оЁеге4. 
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